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Abstract 

The decay chain b B -> D**°£-uX, D**° -» D*+vr-, D*+ -> D°vr+, D° -» (KvrorK37r) is 
identified in a sample of 3.9 million hadronic Z decays collected with the OPAL detector at LEP. 
The branching ratio BR (b -> B) x BR (B -> D? £~PX) x BR (D? -> D*+vr-) is measured to 
be (2.64 ± 0.79 (stat) ± 0.39 (syst)) x 10" 3 for the J p = 1+ (D?) state. For decays into the 
J p = 2 + (T>2°) state, an upper limit of 1.4 x 10 -3 is placed on the branching ratio at the 95% 
confidence level. 
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1 Introduction 



Semileptonic B decays to orbitally-excited P-wave charm mesons (D**) are of interest for several 
reasons. Firstly, they constitute a significant fraction of B semileptonic decays, thereby account- 
ing for some of the difference between the inclusive measurements and the sum of the exclusive B 
decay modes to D*£u and Dlu [3l|T1|2]. They also contribute the major source of systematic error 
in the |V c b| measurement at LEP, as a background to the direct decay B° — ► D* + £~za Finally, 
the measured decay properties can be compared with theoretical HQET predictions [U^. 

The D** mesons (sometimes denoted Dj) are composed of a charm quark and a light quark 
in a state of orbital angular momentum L = 1. In the infinite heavy-quark (charm) mass 
limit the D** system is equivalent to a fixed force center one-body problem. Hence the total 
(spin+orbital) angular momentum of the light quark degrees of freedom (labelled by J q = | or 
|) and the spin of the heavy quark are taken as separately conserved The </<j=§ states 

combine with the heavy quark spin to form two states with J p = 1 + (Di) and J p = 2 + (Drj). In 
the infinite charm mass limit, they can only undergo D-wave decay, and therefore have narrow 
widths For J q =\ we have two states with J p = + and J p = 1 + . The J q =\ states can 
decay via S-wave and are expected to be broad, but their experimental observation is still not 
established So in total, for P-wave mesons, four charged and four neutral D** states are 
predicted. Table ^ summarises the properties of the neutral states. 

In this paper we present a new measurement of semileptonic B decays into the narrow neutral 
V>\ and Dg states. We reconstruct events compatible with the decay chain 

b -> B -» D**°fvX, D**° -» D*+ir-, D*+ -» D tt+, D° -> (K7rorK37r). (1) 

Here, and throughout this paper, B refers to the B° and B~ mesons, and charge conjugate modes 
and reactions are always implied. Also, £ refers to both electrons and muons, and the terms 
"kaon" and "pion" denote the charged particles. 

These results update a previous OPAL analysis 0. Similar measurements have been per- 
formed also by ALEPH 0, DELPHI 0, CLEO and ARGUS [TU]. We note the complemen- 
tarity between the LEP and T(4S) analyses, as in the latter, theoretical uncertainties in the B 
semileptonic decay form factors contribute significantly to the overall systematic error. This is 
less important at LEP, due to the high boost of the B hadron and its decay products. 

2 Detector and Monte Carlo Samples 

The OPAL detector is fully described elsewhere A brief description of the main components 
relevant for this analysis follows. Tracking of charged particles is performed by a central detector, 
consisting of a silicon microvertex detector, a vertex chamber, a jet chamber and z-chambers 1 . 
The central detector is inside a solenoid, which provides a uniform axial magnetic field of 0.435 T. 
The silicon microvertex detector consists of two layers of silicon strip detectors; for most of the 
data used in this paper, the inner layer covered a polar angle range of | cos 6*| < 0.83 and the 
outer layer covered |cos#| < 0.77, with an extended coverage for the data taken after 1996. 
This detector provided only ^-coordinate information for data taken in 1991-1992, and also 
z-coordinate information thereafter. The vertex chamber is a precision drift chamber which 
covers the range |cos#| < 0.95. The jet chamber is a large- volume drift chamber, 4.0 m long 

1 A right handed coordinate system is used, with positive z along the e~ beam direction and x pointing towards 
the center of the LEP ring. The polar and azimuthal angles are denoted by 9 and <f), and the origin is taken to 
be the center of the detector. 
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and 3.7 m in diameter, providing both tracking and ionisation energy loss (dE/dx) information. 
The z-chambers measure the z-coordinate of tracks as they leave the jet chamber in the range 
| cos 9 1 < 0.72. Immediately outside the tracking volume is the solenoid and a time-of- flight 
counter array followed by an electromagnetic shower presampler and a lead-glass electromagnetic 
calorimeter. The return yoke of the magnet lies outside the electromagnetic calorimeter and is 
instrumented with limited streamer chambers and thin gap chambers. It is used as a hadron 
calorimeter and assists in the reconstruction of muons. The outermost part of the detector is 
made up by layers of muon chambers. 

A Monte Carlo simulation of hadronic Z decays of about five times the size of the recorded 
data sample for b flavoured events, and about the same size as the recorded data for the other 
flavours, is used in the analysis. In addition, signal Monte Carlo samples were generated with 
approximately one thousand times the expected number of signal events in the data. The 
simulated events were generated using JETSET 7.4 ^2] with the b and c quark fragmentation 
modelled according to the parameterisation of Peterson et al. A global fit to OPAL data 
has been performed to optimise the JETSET parameters [Tl|- These events were processed using 
the full OPAL detector simulation ^2] and analysed in the same manner as the data. 

3 Event Selection 

This analysis identifies decays compatible with the decay chain of Eq. ^ illustrated in Fig. ^ We 
attempt to identify the 5 relevant tracks in the case of the D° — ► K-7T mode, or 7 tracks for the 
D° — > K 3ir mode, that match the correct event topology. To distinguish between the different 
charged pions involved, we denote as it** the one coming from the D**° decay. Similarly, 7r s i ow 
comes from the D* + meson, its name arising from the small mass difference (Am*) between the 
D* + and the D°, leading to a very low pion momentum in the D* + rest frame. Event selection 
criteria are applied to the data and Monte Carlo samples in five stages, namely hadronic Z 
preselection, reconstruction of D* + ^~candidates, identification of the best 7r**, selection of the 
best overall candidate, and final background suppression cuts specific to the D° — ► K-zr and 
D° — ► K 3tt modes. 

3.1 Hadronic Z Decay Preselection 

Hadronic Z decays collected with the OPAL detector at e + e~ centre-of-mass energies near the Z 
resonance are selected using a standard OPAL hadronic event selection |16j . with an additional 
requirement of at least 7 tracks per event. All the tracks must pass a set of standard quality 
cuts |17j . The selection efficiency for hadronic Z decays is (98.1 ± 0.5)% ^Zj with a background 
of (0.11 ±0.03)%. Only data that were taken with the silicon microvertex detector in operation 
are used in this analysis. After this preselection, the resulting data sample collected in the years 
1991-2000 consists of 3 904 417 events selected from a total integrated luminosity of 180.8 pb _1 
collected at the Z resonance. About 11% of the data used were recorded in the years 1996-2000. 

3.2 Reconstruction of D* + £~ Candidates 

Candidate events must have two or more jets, defined using a cone algorithm |18j . with radius 
of 0.7 radians and 10 GeV minimum energy. The primary vertex of the event is reconstructed 
using a beam spot constraint ^H]- We then identify muon or electron candidates with momentum 
greater than 3 or 2 GeV, respectively. The muon cut is more stringent to reduce the background 
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from pions misidentified as muons. An artificial neural network is used to identify electrons 20 , 
and photon conversions are rejected as described in Muons are selected as described in |17j . 

Next, we look for tracks consistent with the D* + hypothesis and within the same jet as 
the lepton. Firstly, a kaon candidate with the same charge sign as the lepton, and momentum 
greater than 1.0 GeV, is required to have a dE/dx probability greater than 1% for the kaon 
hypothesis. This requirement is increased to 10% if the D* + energy is smaller than half the 
beam energy, where K/V separation is more powerful. Also, the measured dE/dx of the kaon 
candidate is required to be smaller than the expected dE/dx for the pion hypothesis. 

A pion of opposite charge is then sought in order to form a D° — ► K-zr candidate; for D° — > K 3ir, 
two additional (pion) tracks of opposite charge are identified. The reconstructed D° mass peak 
has a resolution of about 25 MeV; a loose cut of 1.79 < M(D°) < 1.94 GeV is applied. A 7r s i ow 
track candidate of opposite charge to the kaon is combined with the D° to form a D* + . The 
mass difference Am* = M(D* + ) — M(D°) forms a sharp peak of about 1 MeV resolution, and 
we require 140.5 < Am* < 149.5 MeV at this stage. The energy of the D* + is required to be 
greater than 15% of the beam energy. A D° vertex is required to be successfully reconstructed 
using the 2 (or 4) candidate tracks. A kinematic fit performed on the D* + tracks using the 
known D* + and D° masses as constraints is required to converge. We note that at this point, 
several different D* + £~ candidates per event can exist. 

3.3 Selection of Best Overall Candidate 

We now proceed to select the best ir** candidate for a given ~D* + £~ assignment, where tt** 
candidates are required at this stage to have an energy greater than 0.6 GeV and the same 
charge as the lepton. A B vertex is constructed using the tracks of the tt**, 7r s i ow , lepton, and 
reconstructed D°. The primary vertex is then recomputed, excluding all of the above 5 (or 7) 
tracks. The new x 2 of the primary vertex is obtained, and for each T)* + £~ combination the 7r** 
candidate that gives the smallest combined x 2 °f the primary and B vertices is chosen. The B 
and D decay lengths are recomputed with the new vertex positions. In order to select only well 
reconstructed B vertices, we require an estimated uncertainty of less than 1 mm on the B decay 
length. 

Having selected the best tt** candidate for the given T)* + £~ assignment, we reconstruct the 
mass difference Am** = M(D**°) - M(D*+), and restrict it to the range 0.14-1.10 GeV. We 
then select for further analysis only one D* + £~ candidate per event, namely the one with highest 
D* + kinematic fit probability. This cut greatly reduces the combinatorial background arising 
from fake combinations when forming the D* + £~ . 

3.4 Final Cuts 

The main backgrounds at this point arise from fragmentation, where a pion from the b quark 
fragmentation fakes a tt** from a D**°, and combinatorial backgrounds; we include in the 
fragmentation category pions from B** — > B7rA. The cuts for reducing these backgrounds 
differ somewhat for the two D° decay modes, the D° — > K 3tt mode having larger backgrounds, 
and correspondingly stronger cuts. 

To reduce the fragmentation background, we make the following requirements: 
• The tt** energy must be greater than 1 GeV. 
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• The output of a neural network applied to the ir** track is required to be greater than 0.6 
for the D° — > K 3-7T mode, and greater than 0.3 for the D° — > Kir mode. This neural network 
uses momentum, transverse momentum to the jet, and impact parameter significance with 
respect to the primary vertex; it is used to distinguish between tracks from the primary 
vertex and genuine B decay tracks [22]. 

• The B decay length significance Ib/&i b is required to be greater than 1.5, where Ib is the 
signed B decay length and ai B its error. It is signed positive if the B vertex is displaced from 
the primary vertex in the same direction as the jet momentum, and negative otherwise. 

• Similarly, Zb-d/o7 B -d > —2.0 for the signed decay length significance between the B and 
D vertices. It is signed positive if the D vertex is displaced from the B vertex in the same 
direction as the jet momentum, and negative otherwise. 

The combinatorial background is suppressed using the following cuts: 

• 143 < Am* < 148 MeV for D° -» K3vr, 142 < Am* < 149 MeV for D° -» Ktt. 

• 1.815 < M(D°) < 1.915 GeV for D° -» K3vr. 

• cos 9* > —0.9, where 9* is the angle between the kaon and the D° boost direction, calculated 
in the D° rest frame. This uses the fact that the D° is a pseudoscalar, whereas the 
background, particularly that resulting from particle misidentification, tends to peak at 
large negative values of cos 9*. 

• For the D° — > K 3tt mode, the product of the neural network applied to the 4 tracks from 
the D° decay, is required to be greater than 0.2. 

Fake leptons and non-B semileptonic decays appear at a smaller rate than the previous 
backgrounds, but their invariant mass distributions can peak in the expected signal region. The 
non-B semileptonic background arises primarily from charm, tau and b-baryon decays, and pion 
misidentification. Therefore we impose the following additional cuts: 

• The D* + £ _ invariant mass is restricted to the range 2.8-5.0 GeV. 

• The lepton momentum transverse to the jet: px t > 0-9 GeV for D° — > K 3n, and > 0.6 GeV 
for D° — > Ktt. This cut is effective in reducing misidentified lepton backgrounds and does 
not reduce the signal significance. 

These selection criteria were designed to maximise the efficiency y. purity using only quantities 
well modelled by the Monte Carlo simulation. As an example, Fig. [2] presents the data versus 
Monte Carlo distributions after all cuts for four of the most relevant quantities used in the 
analysis. Since the final analysis consists of searches for peaks in the distribution of Am**, the 
cuts were also tuned so that the simulated background Am** distribution shows no peak in the 
signal region, 0.35 < Am** < 0.55 GeV. 

4 Backgrounds and Wrong Sign Sample 

Applying the previous selection to our combined D° — > {Ktt or K 3ir) Monte Carlo samples, the 
major background process is B° — > D* + l~v, where the D* + is combined with an additional pion 
of the right sign. This comprises 65% of the total background. The decays B — > D**°£~i?X 
and B — > D** + £~PX, where some of the tracks are incorrectly matched, amount to 16% (using 
the branching fractions measured in this paper). These three distinct physical processes have in 
common that the 7r** track is incorrectly identified, and therefore they constitute the major part 
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of the above mentioned fragmentation background. The kinematics of this dominant background 
should be equivalent to the one obtained by requiring a tt** of the opposite charge from that 
expected for D**° decay. The Monte Carlo Am** distributions for background right sign and 
wrong sign tt** are shown in Fig. |2K,b. They are fitted to a functional form x' y e~^ x where 
x = Am** — m^; the fitted values of the 7 and (5 parameters for the two samples are consistent. 
Fig. |3t shows the OPAL data wrong sign Am** distribution; the fitted 7 and (5 parameters agree 
with those obtained in Monte Carlo. 

Fake leptons, B s , and non-B semileptonic decays constitute less severe backgrounds amount- 
ing to 3%, 5% and 11% of the total background, respectively. We find the contribution from Z 
decays into charm and lighter quarks to be negligible. 

5 Signal Fitting Procedure 

An unbinned maximum likelihood fit is performed simultaneously on the Am** right sign and 
wrong sign distributions. For the right sign fit, the function is given by the sum of two Breit- 
Wigner distributions, each convolved with a Gaussian resolution function, plus the background 
function x^e~^ x described previously; for the wrong sign fit, only the background function is 
used. There are six fitted parameters: the number of and D*, , the normalisation of the right 
sign and wrong sign backgrounds, and the background shape parameters 7 and [5. Thus, the 
background shape and normalisation are obtained directly from data. We fix the narrow states 
masses and widths to the world averages [2], shown in Tabled The detector resolution on the 
reconstructed Am** is fixed at a =8 MeV, as obtained from Monte Carlo. 

To check for biases and to assess the systematic uncertainties, the fitting procedure was 
tested by comparing the results of five simulated experiments, each one with the same statistics 
as the data (so-called "ensemble tests"). The pull distributions for the fitted number of D^, D?, , 
and background events were found to be consistent with zero mean and unit variance for both 
the D° decay modes taken separately as well as when fitting to the two modes combined. The 
uncertainty is taken to be the error on the mean of the pull distributions. We note that the two 
narrow peaks are resolvable and we are able to correctly fit the number of and D?, signal 
events in the Monte Carlo simulations for various input signal branching ratios. 

6 Results 

The final Am** distributions for the D° — > Kir and D° — > K 3tt modes are shown in Fig. ^ 
The numerical results of the fit are also shown, including the 1 a uncertainties on the fitted 
parameters. A signal for the expected narrow V)\ peak is present in each mode separately, with 
a weaker significance for the D° — > Kir mode (8.0±5.0 fitted Dj 1 events) than for D° — > K3tt 
(21.4±7.2 fitted events). No evidence of a D?, signal is present. 

The efficiencies estimated for the state using dedicated signal simulations are found to 
be 7.7 ± 0.4% and 2.3 ± 0.1% for the D° — > Ktt and D° — > K37T modes, respectively, where the 
errors are statistical only. For the D|° state these are 9.2 ± 0.8% and 2.1 ± 0.2%. The product 
branching ratio can then be obtained as: 

BR (b -> B) x BR (B -> D**° I'vX) x BR (D**° -> D* + 7r~) = 

AWD**°)/e D »o m 

(N z /e z ) x R h x 2 x 2 x BR (D*+ -> D vr+) x BR (D° -> Kvr, K 3vr) U 
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The two factors of 2 arise from the two b hadrons in the event, and the two flavors of tagged 
lepton. The fraction, R^, of Z — > bb events in Z hadronic decays, and the branching ratios for 
D*+ -» D°7r+ and D° -» Kvr are taken from 3 ; iV z = 3904417 and e z = (98.1 ± 0.5)%, as 
mentioned in section IXTl 

For decays into the state, we obtain for the D° — ► Ktt mode: 
BR (b -> B) x BR (B -» D? r PX) x BR (D? -» D* + 7r~) = (1.17 ± 0.73 ± 0.27) x 10~ 3 
and for D° -> K 3vr: 

BR (b -> B) x BR (B -» D? r PX) x BR (D? -» D* + 7r~) = (5.30 ± 1.79 ± 0.95) x 10~ 3 

where the first error is statistical and the second systematic. The evaluation of the systematic 
uncertainty is discussed in section [Tj The two product branching ratios measured for the 
D° — > K-7T and D° — > K3tt decay modes agree at the level of 1.9 a. 

The two modes were combined by merging the two data samples and performing a single 
fit, which has the advantage of reducing the systematic error on the background. The resulting 
distribution is shown in Fig[5J The combined efficiency, multiplied by the respective D° decay 
branching ratios, is 0.47% (7.7% x 3.80% + 2.3% x 7.46%). There are 28.7 ± 8.6 events in the 

peak, from which the product branching ratio follows: 

BR (b -» B) x BR (B -» D%fi/X) x BR (D? -> D* + 7r-) = (2.64 ± 0.79 ± 0.39) x 10~ 3 . 

The combined-sample fit yields 3.1 ± 6.9 events for the D*. . The measured branching ratio 
is then consistent with zero: 

BR (b -> B) x BR (B -» D*. f DX) x BR (D*° -> D* + 7r~) = (0.26 ± 0.59 ± 0.35) x 10~ 3 . 

Using the method of Feldman and Cousins this can be converted into a 95% confidence 

level upper limit: 

BR (b -> B) x BR (B D* 2 ° CvX) x BR (D^° -> D* + tt") < 1.39 x 10~ 3 . 

These results were checked in several ways. Firstly, if instead of fixing the V>\ mass we include 
it as an extra parameter to be fitted we obtain: 

BR (b -» B) x BR (B -» D? l~vX) x BR (D? D*+vr") = (2.66 ± 0.80 ± 0.39) x 10" 3 , 

in agreement with the previous result of (2.64±0.79±0.39) x 10" 3 . The fitted mass is (2418. 8± 
5.0) MeV, consistent with the world average of (2422.2 ± 1.9) MeV (Hj. The cuts described in 
section ETH were varied by significant amounts (10%-100%), and the resulting variations in the 
measured branching ratios were found to be compatible with the expected statistical variations. 
To check the procedure for combining the results of the two D° decay modes, we calculated a 
simple weighted average of the two results. This average is consistent with the combined result. 

We also considered the angular decay distribution of the narrow states. In the heavy-quark 
limit, the pure D-wave decay of the J p = 2 + D^ state should be distributed according to 
| sin 2 q, whereas for the J p = 1 + we expect |(1 + 3 cos 2 a) [Jj- This has been observed 
experimentally by CLEO |24j . Here a is the angle between the tt** and the 7r s i ow from the 
D* + decay, evaluated in the rest frame of the D* + . For D^ 1 events we expect more events at 
larger values of | cos a\, and this is evident in Fig. |HK which shows this angle for data and Monte 
Carlo in the signal mass region 0.35 < Am** < 0.55 GeV. The shape of the data distribution 
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is observed to agree with that expected from simulated events, where the signal rate is fixed to 
the one measured here. Fig. m> shows that background events, selected from right sign data in 
the sidebands 0.14 < Am** < 0.3 or 0.6 < Am** < 1.1 GeV, exhibit a flat cos a distribution. 
An enhancement of the resonance peak is expected for higher values of | cosa|, and this is 
confirmed in Fig. Efc as the Am** peak is enhanced with a selection |cosa| > 0.5. Note 
how the overall level of the background drops relative to Fig. [21 Conversely, for | cosa| < 0.5 
(Fig. Ell), the D]* peak greatly diminishes in significance. 

For a D*. enhancement selection, namely | cosa| < 0.2, we find no distinct sharp peak at the 
expected D*, position. This provides additional evidence for the presence of the state and 
absence of the D*, . 

7 Systematic Uncertainties 

The systematic uncertainties in the product branching ratios are shown in Table [21 The 
dominant contributions arise from uncertainties in the background function and the and 
D*, fit parameters. The components of systematic uncertainty shown in the table are: 

• Background function: The dominant uncertainty on the background arises from the pos- 
sibility that some background has a peak in the signal region, therefore biasing the fitted 
number of and D*, signal events. This uncertainty is estimated from the average 
number of (expected-fitted) signal events, as obtained in the Monte Carlo ensemble tests 
described in Section [SI There is also a much smaller background uncertainty on the overall 
background shape. This component of the uncertainty is estimated by comparing the 7, (3 
parameters obtained from Monte Carlo with those obtained from the wrong sign tt** data 
sample. 

• D**° fit parameters: We vary the Dj 1 and D*, masses and widths within the current 
uncertainties [31, and refit the Am** distribution. We also vary the Gaussian resolution 
by ±1 MeV and redo the fit. The resulting variations are added in quadrature. 

All the remaining sources of systematic uncertainty affect only the signal efficiency, as the 
background shape and normalisation are obtained directly from the data. 

• The limited Monte Carlo statistics in the signal samples give rise to a systematic uncer- 
tainty on the estimation of efficiencies. 

• The lepton identification efficiency has an uncertainty of 3% for muons and 4% for electrons 

• The mean and sigma of the normalised Monte Carlo dE/dx distributions were varied by 

±10% mi- 

• Tracking resolution: The systematic uncertainty was assessed in Monte Carlo by applying 
a global 10% degradation to the resolution of all measured track parameters. 

• The lifetime of the B mesons was varied within their measured uncertainty [3]. 

• Theoretical uncertainty in the B — > D**°£~i>X form factors: Different theoretical models 
predict different form factors for the B semileptonic decay, and therefore different ir** 
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energy spectra. This source of systematic uncertainty is dominant at the T(4S) experi- 
ments (§]. Due to the high boost provided in Z — > bb decays, this uncertainty is expected 
to be significantly smaller in the analyses at LEP. This was confirmed by reweighting our 
JETSET ^21 Monte Carlo samples to the forms factors described in |27| . The range of 
parameters describing the form factor calculation was varied within the values constrained 
as in , and also the values needed to approximately reproduce the form factors predicted 
by [221 ■ The maximum variations in the signal efficiency were taken. 

• B fragmentation: The Peterson fragmentation model |13| is used in JETSET 12 to model 
the momentum distribution of the B hadrons. Our Monte Carlo samples were reweighted 
event by event to reproduce the experimental uncertainty on the mean energy of the B 
hadrons. The Peterson parameter was varied within the range obtained in (SHI- The 
variations are in agreement with those obtained using two other fragmentation models |31l 
I32|I33| . again with parameter ranges as determined in 30 . 

• The relevant branching ratios and were varied within the published uncertainties [3]. 

• The hadronic Z preselection efficiency was varied by its uncertainty of ±0.5% |17j . 

8 Conclusions 

We have analysed semileptonic B decays into the narrow Dj 1 and D^ states in events compatible 
with the decay chain b -> B -> B**°£~i?X, D**° -» D*+vr-, D*+ -► D°vr + , D° -> (KvrorK37r). 
The product branching ratio for decays into the state is measured to be: 

BR (b -> B) x BR (B -> D^fPX) x BR (D? -> D* + tt") = (2.64 ± 0.79 ± 0.39) x 10~ 3 , 

where the first error is statistical and the second systematic. We find no evidence for decays 
into the J p = 2 + state, and set a limit on the product branching ratio: 

BR (b -> B) x BR (B -> D* 2 ° fvX) x BR (D* 2 ° -> D* + vr-) < 1.4 x 10" 3 , 

at the 95% confidence level. These results update a previous OPAL analysis [7], and agree with 
similar measurements performed at LEP 0IH1, CLEO [HJ, and ARGUS [TUj . 

Acknowledgements 

We particularly wish to thank the SL Division for the efficient operation of the LEP accelerator 

at all energies and for their close cooperation with our experimental group. In addition to the 

support staff at our own institutions we are pleased to acknowledge 

Support from the scholarship PRAXIS/BD/11460/97, Portugal, 

Department of Energy, USA, 

National Science Foundation, USA, 

Particle Physics and Astronomy Research Council, UK, 

Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and Humanities, 
Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and a grant 
under the MEXT International Science Research Program, 



11 



Japanese Society for the Promotion of Science (JSPS), 
German Israeli Bi-national Science Foundation (GIF), 
Bundesministerium fur Bildung und Forschung, Germany, 
National Research Council of Canada, 

Hungarian Foundation for Scientific Research, OTKA T-029328, and T-038240, 
The NWO/NATO Fund for Scientific Reasearch, the Netherlands. 

References 

[1] J. Richman and P. Burchat, Rev. Mod. Phys. 67, 893 (1995). 

[2] ALEPH Collaboration, D. Buskulic et al, Z. Phys. C73 (1997) 601. 

[3] The Review of Particle Physics, K. Hagiwara et al, Phys. Rev. D66 (2002). 

[4] M. Neubert, Phys. Reports 245, 259 (1994). 

[5] J. Rosner, Comm. Nucl. Part. Phys. 62 (1986) 109. 

[6] S. Godfrey and R. Kokoski, Phys. Rev. D43 (1991) 1679. 

[7] OPAL Collaboration, R. Akers et al, Z. Phys. C67 (1995) 57. 

[8] DELPHI Collaboration, P. Abreu et al, Phys. Lett. B457 (2000) 407. 

[9] CLEO Collaboration, A. Anastassov et al, Phys. Rev. Lett. 80 (1998) 4127. 

[10] ARGUS Collaboration, H. Albrecht et al, Z. Phys. C57 (1993) 533. 

[11] OPAL Collaboration, K. Ahmet et al, Nucl. Instr. Meth. A305 (1991) 275; 
P.P. Allport et al, Nucl. Instr. Meth. A324 (1993) 34; 
P.P. Allport et al, Nucl. Instr. Meth. A346 (1994) 476; 
O. Biebel et al, Nucl. Instr. Meth. A323 (1992) 169; 
M. Hauschild et al, Nucl. Instr. Meth. A314 (1992) 74; 
S. Anderson et al, Nucl. Instr. Meth. A403 (1998) 326. 

[12] T. SjSstrand, Comp. Phys. Comm. 82 (1994) 74. 

[13] C. Peterson, D. Schlatter, I. Schmitt and P. Zerwas, Phys. Rev. D27 (1983) 

[14] OPAL Collaboration, G. Alexander et al, Z. Phys. C69 (1996) 543. 

[15] J. Allison et al, Nucl. Instr. Meth. A305 (1992) 47. 

[16] OPAL Collaboration, G. Alexander et al, Z. Phys. C52 (1991) 175. 

[17] OPAL Collaboration, R. Akers et al, Z. Phys. C65 (1995) 17. 

[18] OPAL Collaboration, R. Akers et al, Z. Phys. C 63 (1994) 197. 

[19] OPAL Collaboration, K. Ackerstaff et al, Z. Phys. C 74 (1997) 1. 

[20] OPAL Collaboration, G. Alexander et al, Z. Phys. C 70 (1996) 357. 

[21] OPAL Collaboration, P. D. Acton et al, Z. Phys. C 58 (1993) 523. 



12 



[22] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C23 (2002) 437. 

[23] G. Feldman and R. Cousins, Phys. Rev. D57 (1998) 3873. 

[24] CLEO Collaboration, P. Avery et al, Phys. Lett. B331 (1994) 236. 

[25] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C8 (1999) 217. 

[26] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C21 (2001) 23. 

[27] A. Leibovich, Z. Ligeti, I. Stewart and M. Wise, Phys. Rev. D57 (1998) 308. 

[28] OPAL Collaboration, G. Abbiendi et al, Phys. Lett. B482 (2000) 15. 

[29] D. Scora and N. Isgur, Phys. Rev. D52 (1995) 2783. 

[30] OPAL Collaboration, G. Abbiendi et al, CERN-EP/2002-051, submitted to Eur. Phys. J. 
C. 

[31] V.G. Kartvelishvili, A.K. Likhoded, V.A. Petrov, Phys. Lett. B78 (1978) 615. 

[32] P.D.B. Collins and T.P Spiller, J. Phys. Gil (1985) 1289. 

[33] M.G. Bowler, Z. Phys. Cll (1981) 169. 



13 









Mass 


Width 


Decay Modes 




J p 




(MeV) 


(MeV) 






1+ 


1/2 


~2470 


> 250 


D*7T 




0+ 


1/2 


~2400 


> 170 


Dvr 


D?(2420) 


1+ 


3/2 


2422.2 ± 1.8 


l8-9±f;I 


D*7T 


Df (2460) 


2+ 


3/2 


2458.9 ± 2.0 


23 ± 5 


Dvr, D*tt 



Table 1: Properties of neutral orbitally-excited charm mesons. The quantum number J q is 
the total angular momentum of the light constituents of the meson. Masses and widths of the 
narrow J q = 3/2 states are experimentally determined jSj- For the light J q = 1/2 states, we 
quote theoretical estimates of the masses and widths [§] • 





Systematic Uncertainty (xlO 3 ) 


Source 


JL v /t muuv 


^Kjr rnonp 

J. \. KJ 1 1 lllv vlV 


(!oiTi ni n pH 

VyvJlllUlllCU. 


Background function 


± 0.23 


± 0.54 


± 0.17 


I l-i 1 V i f"l+" T^P T"51 TTI f^1"PT*Q 
1 2 ^ jJdl CXlllC. LCI O 


+0.09 
-0.06 


+0.54 
-0.42 


+0.22 
-0.17 


Signal simulation statistics 


± 0.06 


± 0.20 


± 0.12 


Lepton ID 


± 0.04 


± 0.18 


± 0.09 


dE/dx 


± 0.05 


± 0.26 


± 0.11 


Tracking resolution 


± 0.04 


± 0.27 


± 0.10 


B lifetime 


± 0.01 


± 0.05 


± 0.03 


Theoretical uncertainty 


± 0.03 


± 0.16 


± 0.08 


B fragmentation 


± 0.04 


± 0.26 


± 0.10 


Rb 


± 0.004 


± 0.02 


± 0.009 


BR(D*+ -» D°7r) 


± 0.008 


± 0.04 


± 0.02 


BR(D° -> Kvr) 


± 0.03 




± 0.02 


BR(D° -> K3vr) 




± 0.22 


± 0.13 


Hadronic Z preselection 


± 0.006 


± 0.03 


± 0.01 


Total 


+0.27 
-0.26 


+0.98 
-0.92 


+0.40 
-0.38 



Table 2: Systematic uncertainties on the product branching ratios. The combined column is 
for the combination resulting from merging the two data samples and performing a single fit. 
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Figure 2: Distributions of data, simulated signal, and backgrounds after all cuts for: a) 
probability of kinematic fit, b) lepton momentum, c) n** energy, d) ir** neural network output. 
Both D° decay modes were combined, and the signal rate is the one measured in this paper. 
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Figure 3: Am** distributions for the combined modes D° — ► (K7rorK37r): a) Monte Carlo 
background right sign events, b) Monte Carlo background wrong sign events, c) OPAL wrong 
sign Am** data. The superimposed lines are the fits to the functional form x^e~^ x with x = 
Am** — m w . 



17 



DjD 2 » 

11 



OPAL 



- 1 — i — i — i — | — i — i — i — r— i — i — i — | — i — i — i — | — i — i — r 

• OPAL D%Ktt 
Dj events = 8.0 ± 5.0 
D 2 * events = 0.9 ± 4.5 
Background events = 118.1 ± 12.1 

7=0.35 ±0.15 
p = 3.1 ± 0.6 



14 
12 
10 
8 
6 
4 
2 




J I I 1— J I I I L 



a)- 




0.2 0.3 0.4 



0.9 1 

** 

Am (GeV) 



1.1 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — | — i — i — r 




• OPAL D u ^K3tc 
Dj events = 21.4 ± 7.2 
D 2 * events = 1.6 ± 5.2 
Background events = 96.0 ± 11.4 
y=0.59 ±0.19 
p = 5.0 ± 0.8 



b) 



J. 



0.2 0.3 0.4 0.5 0.6 0.7 



0.8 0.9 1 1.1 

** 

Am (GeV) 



Figure 4: Am** distribution for a) the D° — > Ktt mode, b) the D° — > K 3tt mode. The 
superimposed lines show the overall fit and the background fitted shape. The expected positions 
of the D]* and D*, are indicated by the arrows. 
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Figure 5: Am** distribution for D — > Ktt and D — > K 3ir combined. The superimposed lines 
show the overall fit and the background fitted shape. The expected positions of the and D*. 
are indicated by the arrows. 
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Figure 6: Distributions dependent on the decay angle a: a) | cos(a)| for the combined D° — ► K-zr 
and D° ^K37r data and Monte Carlo in the 0.35 < Am** < 0.55 GeV signal band; b) 
complementary events to a) in the 0.14 < Am** < 0.3 and 0.6 < Am** < 1.1 GeV sidebands; 
c) Am** data distributions for D° — > Ktt and D° — > K3-7T combined data in the | cos(a)| > 0.5 
region (expected to enhance the Di signal); d) complementary events to c) for |cos(q)| < 0.5 
(expected to suppress the Di signal). The expected positions of the Di and D*. are indicated by 
arrows. 
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